Photo control of transport properties in disorderd wire; 
average conductance, conductance statistics, and time-reversal symmetry 



Takuya Kitagawa , Takashi Oka 1 ' 2 , Eu gene Demler 
1 Physics Department, Harvard University, Cambridge, 
Massachusetts 02138, USA, 2 Department of Physics, 
Faculty of Science, University of Tokyo, Tokyo 113-0033, Japan 
(Dated: January 4, 2012) 

In this paper, we study the full conductance statistics of disordered one dimensional wire under 
the application of light. We develop the transfer matrix method for periodically driven systems 
to analyze the conductance of large system with small frequency of light, where coherent photon 
absorptions play important role to determine not only the average but also the shape of conductance 
distributions. The average conductance under the application of light results from the competition 
between dynamic localization and effective dimension increase, and shows non-monotonic behav- 
ior as a function of driving amplitude. On the other hand, the shape of conductance distribution 
displays crossover phenomena in the intermediate disorder strength; the application of light dramat- 
ically changes the distribution from log-normal to normal distributions. Furthermore, we propose 
that conductance of disordered systems can be controlled by engineering the shape, frequency and 
amplitude of light. Change of the shape of driving field controls the time-reversals symmetry and 
the disordered system shows analogous behavior as negative magneto-resistance known in static 
weak localization. A small change of frequency and amplitude of light leads to a large change of 
conductance, displaying giant-opto response. Our work advances the perspective to control the 
mean as well as the full conductance statistics by coherently driving disordered systems. 



I. INTRODUCTION 

The study of disordered systems has been one of the 
central subjects in physics not only because its direct 
relevance to the material science and technology, but 
also because it leads to rich fundamental physics as first 
pointed out by Anderson pQ. While Anderson first con- 
sidered the problem of non-interacting electrons hopping 
on a square lattice with on-site disorders^, the study was 
extended to include the effect of inelastic collisions with 
phonons^El anc [ among particles^, to different geometries 
such as hexagonal lattice in graphene^H and to different 
physical s ystem s such as classical lightP^ and neutral 
cold atomsSSEU, Such extensions not only significantly 
advanced the understanding of disordered systems and 
allowed the observation of Anderson localization, but also 
they brought many surprises, such as the phenomenon of 
anti-weak localization in graphene^ and the presence of 
insulator-metal transitions in interacting systems^. In 
this paper, we further explore the rich physics of disor- 
dered systems by investigating it under the application 
of coherent periodic drive such as light. 

Disordered systems in one dimension have been inten- 
sively studied in the past, and it has been shown through 
numerical as well as analytical meth ods that all the states 
in one dimension are localize d 12 * 13 *. For a finite system, 
the magnitude of average conductance is determined by 
the ratio of localization length £ and system size L; for 
£/L <C 1, the conductance is exponentially suppressed 
and the system is an insulator; on the other hand, for 
%/L ^> 1 the system behaves as a conductor. In this pa- 
per, we demonstrate that the application of light on such 
disordered systems can dramatically change the average 
conductance. 



There are two competing effects that alter the aver- 
age conductance. One is the effect of dimensional in- 
crease. In strongly disordered systems, photons can be 
thought of as providing energies to hop from one localized 
state to another by bridging the energy gaps between the 
localized states. Naively, such photon-assisted hopping 
decreases the effect of disorders and increases the con- 
ductance. We show that this photon-assisted hopping 
provides the effective dimensional increase by inducing 
additional conducting channels called Floquet channels, 
which in return results in the increase of the conductance. 
Previously, using such a picture, a pioneering work by 
Martinez and Molina [14] studied the increase of localiza- 
tion length through the low frequency drive . The sec- 
ond effect is called dynamic localizatio n^ 5 * 16 *, which de- 
creases conductance when disorders are weak and driving 
is strong. This effect has been experimentally observed in 
the context of semiconductor superlattices^ and neutral 
cold atoms^. The origin of dynamic localization is the 
decrease of time-averaged hopping amplitudes for strong 
AC driving of the system ^ * 19 * 20 * . For generic strength 
of disorders, these two effects of dimension increase and 
dynamical localization compete, and as a result, the av- 
erage conductance displays non-monotonic behavior as a 
function of driving amplitude^. 

However, the average conductance is only a partial 
story; in static disordered systems, it is known that 
higher moments displays rich physics such a s the uni- 
versal fluctuations for weakly localized regime^ 2 ^, and 
change of the shape of conductance distributions for in- 
creasing disorder strength. In this paper, we study the 
full conductance distributions to investigate how, for ex- 
ample, the fluctuations of conductance changes upon the 
application of light. Recent work by Gopar and MolinsP^ 
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investigated the behavior of conductance distributions 
for small system size L = 10, and they showed that 
the distribution can be modified through a sharp cut- 
off in the high frequency regime. Their study employed 
so-called continued fra ction m ethod to compute the Flo- 
quet Green's function J 14 * 23 * 24 l This method, being the 
solution that uses the full matrix for the Hilbert space of 
L sites, is restricted to small system size. 

Light application is expected to yield a large effect 
when photon-assisted transport plays an important role, 
where spatial disorders are avoided through the absorp- 
tions and emissions of photons. Such photon-assisted 
transport occurs when two sites in the system are coupled 
resonantly through light. For small system size L = 10 
studied in the previous workP^l, a large driving frequency 
and amplitude compared to the finite size gap AE ~1/L 
were required to satisfy the resonant condition. The mag- 
nitude of driving amplitude suggested by this work is 
challenging to achieve in experiments. 

In this paper, we demonstrate that a large effect of 
photon-assisted transport appears even for small fre- 
quency and small amplitude of driving for large system 
size, resulting in the giant opto-response of the conduc- 
tance in disordered systems. In order to study this in- 
teresting regime, we develop a transfer matrix method 
adapted for periodically driven systems. This method 
allows us to study L = 601 sites when ~ 10 photons 
involve in the transport process. Using this formal- 
ism, we study the conductance distributions for various 
intensity of light and for various strength of disorders 
to give comprehensive picture of conductance behaviors 
across localized and conducting regimes. The analysis 
of conductance distributions in this regime reveal novel 
physics that was not observed before; a dramatic change 
of the shape of conductance distributions in the inter- 
mediate disorder regime, which results in the crossover 
from log-normal like distribution to normal distributions. 
This result opens the perspective to change the conduc- 
tance statistics in disordered systems through the coher- 
ent driving. 

From the point of view of conductance control, ac- 
tive coherent drive through light allows one to effectively 
change parameters of the system. In this paper, we pro- 
pose the possibility of changing the symmetry class of 
driven systems; by choosing the appropriate shape of the 
driving, one can preserve or break the time reversal sym- 
metry. Time-reversal symmetry plays an important role 
to determine the conductance in the static systems, and 
we expect similar effects should exist in the driven sys- 
tems. Here we demonstrate that such control of time- 
reversal symmetry indeed changes the conductance prop- 
erties in the disordered systems. 

Moreover, in experimental and practical situations, it 
is interesting to ask if and how much conductance can be 
controlled by changing the frequency as well as the inten- 
sity of light for a sample with a given realization of disor- 
ders. Because only slightly different frequency and inten- 
sity of light couples very different localized states due to 



the dense density of states of disordered systems, a small 
change of driving frequency can induce a large change in 
the conductance. We will demonstrate, through the for- 
malisms developed in this paper, that the conductance 
of a sample with a given realization of disorders can be 
changed by large amount through the small change of 
frequency and intensity of light. 

While in this paper, we use the language of con- 
densed matter physics, the same consideration can be 
applied to cold atom systems. Recent experiments 
have studied the transport and localization properties 
in one and three dimensional disordered cold atom 
systems ^ * 11 * 25 *. Moreo ver, c oherent AC drive through 
shaking is demonstrated^^. Due to the versatile con- 
trols and probes in cold atom systesms, they provide a 
promising candidate system to confirm the predictions 
made in this paper. 

This paper is organized as follows. In Section [TTJ we 
summarize the main finding of this paper. After a brief 
review of the results of conductance distributions in static 
disordered systems in Section [II A[ we describe the sys- 
tems considered in this paper in Section |IIB| Here we 
also discuss the central concepts of this paper, dynamic 
localizations and dimensional increase, through the Flo- 
quet formalism. In Section II C[ we give the main re- 
sults of numerical calculations for the average (Fig. [2| 
and full distribution (Fig. [3| of conductance, which dis- 
play different behaviors depending on whether the sys- 
tem is in strongly localized regime, intermediate regime 
or conducting regime. The detailed formalism of Floquet 
transfer matrix which allowed the evaluation of conduc- 
tance distributions for large system size is developed in 



Section |III| In Section [TVl we supplement the study pre- 
sented in Fig. [2] and Fig. [3] through the comparison of av- 
erage, fluctuation and skewness of conductance distribu- 
tions for various strength of light intensities and disorder. 
Moreover, we compare these moments of conductance in 
the driven systems with those of static quasi- ID systems 
for various number of channels and strength of disorder. 
In Section [IVB I and Section |IVC[ we analyze how active 
control of light leads to the manipulations of conductance 
in disordered systems. As we discussed above, we present 



in Section [IV B| that, by choosing the appropriate shape 
of the light, one can either break or preserve the time 
reversal symmetry. In Section |IV C[ we study how the 
conductance of a sample with a given realization of dis- 
order depends on the frequency and intensity of light. 
We demonstrate that a large change in the conductance 
can be induced by a small change of either frequency or 
intensity of light. 



II. SUMMARY OF RESULTS 

A. Static disordered systems 

The physics of conductance in the disordered systems 
under the application of light can be best illustrated by 
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contrasting with the conductance in the static limit, or 
the limit in which the intensity of light is set to zero. 
Thus we first summarize the known results in the static 
systems. 

While the study of disordered systems is challenging 
in dimensions higher than one, elegant analysis through 
random matrix theory allows the understanding of not 
only the average value of conducntace, but also their full 
conductance distributions for non-interacting particles in 
quasi-lD disordered systems^. In particular, the solu- 
tion of Dorokhov-Mello-Pereyra-Kumar (DMPK) equa- 
tion gives a reliable, analytic understanding of conduc- 
tance distributions in quasi- ID systems^. Here quasi- 
lD means the system is one dimension but has width W 
which is much smaller than the length of the system L, 
i.e. W <C L. Such analysis shows that a single dimen- 
sionless parameter a = jl/L determines the conductance 
distributions in the two extreme limit of jl/L <C 1 and 
jl/L ^> 1, where I is mean-free path, L is the system 
size and 7 = j3N + 2 — f3 depends on the number of 
channels/the width of the systems N . (3 characterizes 
the symmetry of the systems, and takes = 1 in the 
presence of time-reversal symmetry with no spin-orbit 
coupling, and f3 = 2 in the absence of time-reversal sym- 
metry. When jl/L <C 1, quasi- ID system is in local- 
ized regime, and the average conductance is exponen- 
tially suppressed, i.e. (G/Gq) oc exp (— L/(2jl)) as a 
function of system size L. Here Go = 2e 2 jh is the natu- 
ral unit of conductance. In this expression, £ = 7Z can be 
seen to correspond to localization length. Moreover, the 
distribution of the conductance becomes log-normal, or, 
the distribution of log of the conductance — ln(G/Go) be- 
comes Gaussian with mean value — (ln(G/Go)) = 2L/jl 
and its fluctuations Var [ln(G/Go)] = 4L/jl such that 

P(HG/G )) oc exp|-(-ln(G/G )-^) 2 /(f)}- 

In the other limit jl/L ^> 1, the system shows Drude- 
type conducting behaviors with average conductance 
(G/Go) = Nl/L. The distributions becomes Gaussian 
and the fluctuation of G/Gp displays the universal fea- 

Jm n 



tures Var [G/G ] 
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Therefore, conductances of 



quasi- ID systems increase as one increases the width 
TV, and in particular, for sufficiently large TV, there is 
a crossover from a insulating regime to a conducting 
regime. In analogy with this static results, we show be- 
low that the application of light effectively increases the 
number of channels N in the system, which lead to sim- 
ilar crossover behaviors. 

We note that the localization length in the strongly 
disordered regime crucially depends on the symmetry of 
the system such as time-reversal symmetry, which deter- 
mines the value of f3. For N > 1, the localization length 
becomes longer when the time-reversal symmetry is bro- 
ken. Moreover, the weak localization in higher dimen- 
sions is known to display negative magneto-resistance 
since magnetic field breaks time-reversal symmetry and 
this destroys the constructive interference between time- 
reversal paths. Thus, time-reversal symmetry controls 
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FIG. 1: a) Illustration of the disordered one dimensional sys- 
tems under the application of light considered in this paper. 
Sites that are driven by light is colored as red, and sites that 
are not driven by light act as leads and colored white. Dis- 
orders are also present at sites colored as red, but absent in 
the white regions corresponding to the leads, b) The one di- 
mensional periodically driven system as illustrated in a) can 
be understood in terms of two dimensional time-independent 
Hamiltonian of Eq. Q with additional dimension given by 
energy ladders. A photon can be absorbed or emittted in the 
system, leading to the hopping in the energy ladders. We 
illustrate in the figure only zero, one and two photon absorp- 
tion processes, but as the intensity of light increases, there are 
processes which absorbs or emits more photons. The linear 
potential in energy direction coming from the energy differ- 
ence of different photon number sectors is illustrated by the 
color gradient in the figure. 



the transport behaviors in disordered systems. We will 
see in Section [IV B| that time-reversal symmetry can be 
controlled through the shape of the driving fields in op- 
tically driven systems. 



B. Systems 

In the following, we consider a tight-binding model for 
spin-less, non-interacting fermions in a one dimensional 
lattice with on-site disorder potentials Wi which takes 
random values between — W/2 and W/2. The model can 
be easily extended to include spin degrees of freedom. 
We incorporate the effect of the application of light as 
a phase accumulation of the hopping. See Fig. [T^i for 
illustration. The Hamiltonian is given by 



H(t) 



-J 



<ij> 



c\ Cj e iAi 



(*) 



(1) 



Here Aij(t) = eAasm(Qt)(i— j) where A is the amplitude 
of vector field related to the electric field strength as E = 
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AQ, a is the lattice constant and e is electron charge. In 
the following, we use the units in which a = 1 and e = 
1. Note that for one dimensional lattice with no width, 
polarization of light does not matter. While we assumed 
that the width of the system is 1 and thus the number 
of channels is equal to 1 above, it is straightforward to 
extend the formalism to the systems with more than one 
channel. 

Now we consider the Schrodinger equation for this 
time-dependent Hamiltonian. We take the Fourier 
transform of the state which separates the pho- 
ton absorptions and quasi-energy, i.e. \^(t)) = 
I-Q% dEe ~ it{E+mQ) \^E(rn)). After the Fourier 
transform, the Schrodinger equation becomes time- 
independent; 

(E + mfi)|^(m)) = H ( m ~ k )\^E(k)) (2) 

k 

where H(m) = ^ J^ /Q dte imQt H(t) is mth Floquet 
component of the Hamiltonian. The state ^^(m)) rep- 
resents the state with quasi-energy E which has absorbed 
m number of photons. Eq. Q encapsulates the idea of 
dimensional increase due to the periodic drive, where 
H(m) acts as a hopping to mth neighbor in the en- 
ergy direction. We illustrate this picture in Fig. [TJd. In 
the presence of light with frequency the energy of 
the system is not conserved but it can change by inte- 
ger multiples of Q. This number of absorbed/emitted 
photons m gives the extra index of states and the time 
evolution of states in this one dimensional driven sys- 
tem is effectively the evolution of a particle in two di- 
mensional lattice where one dimension comes from spa- 
tial dimension and another from energy dimension. In 
our system, H(m) = -J m Y,i>j 4 c j ~ J -m J2i<j 4 c j 
where J m = JB m (eAa) with B m (x) is mth Bessel func- 
tion. Here we used the property of Bessel function that 
B m (x) = B- m (-x). 

One important physics in Eq. ([2| is the presence of 
mVt in the left hand side, where energy difference be- 
tween different photon number sectors acts as a linear 
potential in the energy direction. This term leads to the 
Wannier Stark localization and " Bloch oscillation" in the 
energy direction and localizes the particle in a finite num- 
ber of photon sectors. Thus, photon absorptions gener- 
ically lead to increase of conducting channels through 
increasing photon sectors that are involved in the trans- 
port. We refer to these channels as "Floquet channels" 
to distinguish them from channels of static quasi- ID sys- 
tems. The localization length of the Bloch oscillation 
determines how many Floquet channels are activated for 
a given frequency and intensity of light. 

A rough estimate of the number of Floquet channels in- 
volved in the transport can be given in the weakly driven 
systems with A <C 1. In the absence of disorder, Hamil- 
tonian becomes diagonal in quasi-momentum space and 



takes the form 

H = -23 cos(k - A sin(flt)) 4 Ca- 
fe 

For small A, the hopping amplitude Je in the energy 
direction for a given quasi momentum k is approxi- 
mated by Je — — J sm(k)Bi(A). From the study of 
Bloch oscillations, we know the localization length is 
given by Je/& ~ JB\{A)/£l. For the frequency of 
drive Q = 0.01J we consider in this paper, the num- 
ber of Floquet channels involved in the transport for 
the intensity of light A = 0.1 ~ 0.2 is estimated to be 
Nf = JBi(A)/Q ~ 5 ~ 10. From the estimate above, 
we note that the localization length increases as the fre- 
quency of drive decreases and intensity of light increases. 

Another important effect of the application of light is 
the behavior of zero-photon absorption/emission hopping 
Jo in Fig. [TJd. Since Jo oc Bo(eAa), the hopping ampli- 
tude decreases as the intensity of light is incre ased. This 
is the origin of dynamic localizatio n 1 14 * 15 * 1 ^^. The con- 
ductance of the system is determined by the competition 
between the increase of dimension due to the increase in 
the hopping amplitude in the energy directions J m with 
m/0 and the dynamic localization due to the decrease 
in the hopping amplitude Jq. 

The result of such competitions appears in the in- 
crease/decrease of conductances in a finite wire. Here 
we consider the attachments of leads to the periodically 
driven systems and study the DC response current for a 
small chemical potential difference between left and right 
leads, see Fig. [T] We suppose that light and disorder ex- 
ist from site i = 1 to i = L and refer to this part of 
the system as "central system." The left lead goes from 
i = — oo to i = and the right lead from i = L + 1 to 
i = oo. Here, we take a simple hopping Hamiltonian in 
the leads with hopping strength J^; 

H L = -J L c\ Cj (3) 

<ij> 

In the following, we assume the same strength of the hop- 
ping amplitude Jl in both the left and right leads. We 
emphasize that attachment of leads is conceptually cru- 
cial for determining the occupation numbers of fermions 
in Floquet states. In particular, conductances cannot be 
calculated through Kubo's formula often used in static 
systems unless one makes assumptions about the occu- 
pation number of fermions. 

For a periodically driven system, the relation between 
DC current and the transmission matrix i is knowrP^H, 
and given by 

J ~ ^pump H~ ^res (4) 

Jres = ^|t(£m,-E )| 2 (^L-W?) (5) 

m 

DC current can be separated into two parts. One is a 
pump current J pU mp that exists even when the chemical 
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potentials of the left lead hl and right lead fiR are the 
same. When the system possesses inversion symmetry, 
it is necessarily zero. For our disorder system, J pump is 
generally non-zero for each realization of disorders, but it 
is zero on average and we ignore the contribution in this 
paper. Response current J res depends on the transmis- 
sion amplitude t(E m , Eq) of the particle which enters the 
central system from left lead with with energy Eq = jiL 
and leaves to the right lead with energy E m = fi L -\- mQ. 
This is a natural extension of scattering formalism to 
the systems that allow absorptions and emissions of in- 
teger multiples of the photon energy. Here different pho- 
ton number sectors act as different channels of conduc- 
tions. If we represent t(E m ,Eo) as a matrix t in terms 
of the photon numbers, then the response conductance 
G = Jres/Oi, - I^r) is given by G = G {&t) 0Q - This 
restriction to 00 component is the physical consequence 
of the fact that the fermions coming from leads are al- 
ways at the chemical potential hl- Note the important 
difference from the expression of conductances for multi- 
channel static systems which sums over all the channels, 
i.e. g = GoTr (iH). We calculate the transmission am- 
plitudes i(E m ,Eo) through the transfer matrix method 
adapted for periodically driven systems. In the following 
section, we discuss the main results of the analysis and 
leave the formalism to Section Hill 



C. Main results 

In this paper, we study the conductance distributions 
for large system size L = 601 and small driving frequency 
ft = 0.01 J. In order to study better the behaviors of con- 
ducting regimes, we choose the number of (static) chan- 
nels in the central systems to be 2 by considering the 
quasi- ID systems with width = 2. This makes the max- 
imum possible conductance to be 2Go in the absence of 
disorder. The conductance distributions are calculated 
for various disorder strength and intensity of light by 
repeating the conductance calculations for different re- 
alizations of disorders for 10000 times. In the numeri- 
cal evaluation, it is necessary to terminate the Floquet 
channels to finite numbers. We estimated the number of 
Floquet channels involved in the transport to be ~ 10 in 
Section [II B | and here we retained 20 Floquet channels to 
ensure the convergence of conductance. 



1. Average conductance 

We first analyze the behavior of average conductance 
for various intensity of light and disorder strength. Fig. [2] 
shows the color plot of the average conductance (G/Go) 
in log scale. The behavior of average conductance as a 
function of intensity of light, given by A = eAa, strongly 
depends on the disorder strength. 

For weak disorder strength W = 0.2 J, the localiza- 
tion length is longer than the system size L and the 



conductance is large (G/Go) ~ 1. In this regime, ap- 
plication of light leads to decrease of hopping amplitude 
Jo which is proportional to the zeroth Bessel function , 
i.e. Jo oc Bo (A), lea ding t o dynamic localization (see the 
discussion in Section II B). Thus, the conductance tends 



to decrease as A is increased as one can see left bottom 
of Fig. [2] Moreover, the conductance displays oscilla- 
tions as a function of A. Such conductance oscillation 
has been previously studied in the absence of disorder 
under AC drive 32 , and they have shown that the conduc- 
tance becomes proportional to Bo(AL) 2 for large driving 
frequency. Thus, oscillation period is determined by the 
system size L. We note that the effect of the dynamic 
localization is surprisingly large, considering that zero- 
photon hopping is changed only by a small amount at the 
intensity of light considered here i.e. the zeroth Bessel 
function takes Bq(A = 0.2) = 0.99. The dependence 
of conductance on Bo(AL) 2 in the limit of zero disorder 
and high frequency suggests that application of coher- 
ent light can decrease the conductance more effectively 
than the naive guess coming from zero-photon hopping 
strength Jq oc Bq(A). Detailed analytical study of such 
giant dynamic localization in the disordered systems is 
an interesting future work. 



On the other hand, for intermediate disorder of W = 
0.6J, the average conductance increases as A increases 
for small A, see middle bottom of Fig. [2] This increase 
can be understood from the increase of Floquet channels 



as we have explained in Section |HB] However, as the light 
intensity is further increased, the increase of average con- 
ductance stops and starts decreasing around A ~ 0.05. 
Because A is still small, the number of Floquet channels 
is increasing for increasing A, and thus, such decrease 
comes from the dynamic localization, where zero-photon 
absorption/emission hopping amplitude Jo is decreased. 
We note that the change of average conductance is rel- 
atively small, and the average conductance at A = is 
smaller than that of A = 0.18 by about 10%. However 
we will see below that the full conductance distribution 
shows a dramatic change from v 4 = 0to^4 = 0.18. 



For strongly disordered system with W = 1.8 J, the 
change of average conductance due to dimensional in- 
crease, or, the increase of the number of channels could 
be dramatic. In this regime, the conductance G/Go is 
exponentially suppressed, so it is more appropriate to 
consider log (G/Go). In the right bottom of Fig. |2j we 
have plotted the average of log(G/Go) as a function of 
A. We see that (log(G/Go)) increases from w —33 at 
A = to « —15 at A = 0.1. As in the case of inter- 
mediate disorder strength, we see the increase of average 
conductance for small A and decrease for larger A. Such 
non-monotonic behaviors of conductance are suggested 
by previous study of localization length 14 where they cal- 
culated a similar behavior for the localization length as 
a function of amplitude of driving. 
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FIG. 2: (Top) color plot of average conductance (G/Go) for various intensity of light A = eAa and disorder strength W. Here 
the color is assigned in the log scale. The system size is taken to be L = 601 with width (the number of channels) equals 
2, and the driving frequency is Q = 0.01J. (Bottom) The line cut of the average conductance for three different disorder 
strength, W = 0.2, 0.6, 1.8 J, are shown at the bottom, corresponding to weak, intermediate and strongly localized regimes. For 
weakly disordered regime W = 0.2 J (left), the average conductance decreases for increasing intensity of light A due to dynamic 
localization. For intermediate and strongly disordered regimes W = 0.6, 1.8J(middle and right), the average conductance 
increases for small A due to the increase of the number of Floquet channels. For larger A, dynamic localization becomes 
stronger than the increase of Floquet channels, and the average conductance starts decreasing. For W = 1.8 J, we plot the 
average of logarithmic conductance (\og(G / Go)} appropriate for strongly localized regime. 



2. Full conductance distribution 

The study of static disordered systems shows rich 
physics not only in the average conductance but also in 
the full conductance distributions. Motivated by this re- 
sult, we study the full conductance distributions under 
the application of light. In Fig. [3J we have plotted the 
conductance distributions for various strength of disorder 
and intensity of light. 

For weak disorder strength W = 0.2 J and zero inten- 
sity of light (static limit), we see, in agreement with the 
solution of DMPK equation, that the conductance distri- 
bution is normal distribution (see bottom left of Fig. [3]). 
As the intensity of light A is increased from zero, the av- 
erage conductance decreases due to dynamic localization, 
as we have seen above. While the shape of distribution 



stays normal, we see that the width of the distribution, 
or the fluctuation of the conductance distributions be- 
comes much smaller. The square root of the variance 
of G/G at A = is ^Vaxfc/Go] ~ 0.32 whereas at 
A = 0.15 is y/Var[G/G ] « 0.08. Considering that the 
average conductance only changes slightly, which goes 
from (G/G ) « 0.99 at A = to (G/G ) « 0.87 at 
A = 0.15, the application of light has more dramatic 
effects on the fluctuation of the conductance. 

The behavior of average conductance is quite different 
for weakly disordered regimes and strongly disordered 
regime where the former decreases as a function of the 
intensity of light A, the latter increases. Nonetheless, 
strongly disordered regime with W = 1.8 J displays a 
similar behaviors of the decrease of fluctuations as the 
intensity of light A is increased, see the right bottom 
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FIG. 3: Full conductance distributions of disordered ID systems under the application of light for disorder strength W = 
0.2J(left), W = 0.6J(middle), W = 1.8J(right) for various intensity of light. Each colored marker on the diagram indicates 
the disorder strength and intensity of light, and corresponding conductance distribution is plotted with the same color. For 
weakly disordered regime (bottom left), the conductance distributions take normal distributions, and the shape does not change 
upon the application of light. The average conductance decreases, and moreover, the fluctuation significantly decreases as the 
intensity of light A is increased. For strongly disordered regime (bottom right), the distribution of logarithm of conductance is 
plotted. In accord with DMPK equation, the conductance takes log-normal distributions, and this shape does not change upon 
the application of light. As opposed to the weakly disordered regime, the average log-conductance increases as A is increased. 
The most dramatic change occurs in the intermediate regime, where the conductance distribution is strongly skewed akin to 
log-normal distribution in the absence of light. As the light is turned on, the distribution immediately changes its shape and 
approaches normal distribution for larger intensity of light A, displaying crossover behavior. 



of Fig. [3] Note that here we plot the distribution of 
log(G/Go) and not, as in the weakly disordered regime, 
G/Go- In accordance with the expectation from DMPK 
equation, the conductance distribution takes log-normal 
distributions in the static limit. The shape of distribution 
does not change as the intensity of light A is increased. 
While the average of log conductance log G/Go increases, 
the fluctuation of log G/ Go decreases. 

The most dramatic change of conductance distribution 
occurs in the intermediate disorder regime. For zero in- 
tensity of light, the distribution of G/Go shows a distri- 
bution which has a long tail into the large conductance, a 
signature of log- normal like distribution. As the intensity 
of light is increased, the shape of the distribution imme- 
diately changes. It is notable that even for a very small 
intensity of light A = 0.01, the part of the distribution 
with small conductance is strongly suppressed and the 
peak shifts from G/Go to a finite value. As the intensity 
of light is further increased, the distribution approaches 
normal distribution while the fluctuation also decreases. 



We contrast this behavior of full conductance distribu- 
tion with that for the static quasi- ID disorder systems 
with varying number of channels. 

In Fig. [I] a), we have plotted the full conductance dis- 
tribution function of the static quasi- ID systems with 
disorder strength W = 0.4J for varying number of chan- 
nels (width) TV of the system. We see that as the num- 
ber of conducting channels increases, the average con- 
ductance increases by large amount, which goes from 
(G/Go) ~ 0.27 at TV = 1 to (G/G ) ~ 1.32 at TV = 10. 
Such increase of average conductance is accompanied by 
the dramatic change of the shape, where the distribution 
is strongly skewed at TV = 1 but approaches normal dis- 
tribution at TV = 10. Thus as is expected from DMPK 
equation, the system crossovers from insulating regime 
for small number of channels to conducting regimes for 
larger number of channels. 

We juxtaposed this distributions for static systems 
with the conductance distributions under the application 
of light for varying intensity of light in Fig. [4]b). As we 
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FIG. 4: a) Conductance distribution of static quasi-lD dis- 
ordered system for W = 0.4 J and various number of chan- 
nels (width). As the number of channels increases, the av- 
erage conductance increases, and the distribution approaches 
from log-normal to normal distributions. It displays a clear 
crossover behavior from insulating (log-normal) to conduct- 
ing (normal), b) Conductance distribution of driven systems 
with W = 0.6 J for various intensity of light A. While the av- 
erage conductance does not show a large change, the shape of 
the distribution changes as A is increased. This crossover be- 
havior of the shape is quite similar to the crossover observed 
in the static systems for the increasing number of channels. 



have argued in Section |IIB[ the increase of intensity of 
light A increases the number of Floquet channels. This 
picture is analogous to the increase of physical channels 
in the static systems. 

The physical consequence of increase of Floquet chan- 
nels has similarities as well as differences, compared to 
that of increase of physical channels in the static sys- 
tems. First of all, such increase of Floquet channels leads 
to only a moderate increase of average conductance. In 
fact, as we have seen above, the increase of the average 
conductance saturates around A ~ 0.05, and it starts de- 
creasing for larger value of A. As a result, the average 
takes (G/Gq) ~ 0.146 at A = and only slightly larger 
value of (G/Gq) « 0.177 at A = 0.18. On the other 
hand, the change in the shape of distributions is dra- 
matic, and displays the similar crossover behavior. Thus 
the increase of Floquet channels in this regime affects 
the shape of the conductance distribution and induces 
a crossover from log-normal like distribution to normal 
distributions. 

While the change of average conductance was moder- 
ate for the parameters of our study, it is plausible that 
for larger number of physical channels, the light applica- 
tion can induce larger change. Here we considered N = 2 
case, and this implies the maximum conductance possible 
is 2. Increase of the number of channels would increase 
the maximum conductance, and allow larger changes of 
average conductance. 

Our study of conductance distributions under the ap- 
plication of light has a number of extensions. First of all, 
the idea of dimensional increase and increase of Floquet 
channels is general, and applies to any dimension. Static 



disordered systems are known to be localized in 2D but 
disordered systems in 3D are believe to have transitions 
from localized states to delocalized states, depending on 
whether the energy of the state is below or above mo- 
bility edges. Thus, introduction of coherent light on 2D 
disordered system can have larger effects on the average 
conductance, and it is an interesting future work to study 
this physics. 



III. FLOQUET TRANSFER MATRIX METHOD 

In this section, we give the details of transfer matrix 
method adapted for periodically driven systems. As in 
Section |IIB[ we choose the gauge in which light applica- 
tion gives rise to the hopping of a particle with absorp- 
tion/emission of energy. In this gauge, when carrying 
out numerical analysis, one often has to terminate the 
expansion of e lAij ^ in a finite harmonics of driving fre- 
quency Q. Equivalently, one can use the gauge in which 
the light field gives AC electric field E ^ xi cos(fit)c|cj, 
which is computationally more convenient for a finite sys- 
tem because there is only the first harmonics of Q. The 
formalism can be straightforwardly applied to this gauge 
as well, and we describe the details in the Appendix. 

The transmission amplitude t relates the amplitudes of 
incoming modes from the left to those of outgoing modes 
to the right, see Fig. [5] In the leads, there is no disor- 
der or driving fields, so free propagating mode e ±%kx is 
the stationary states with energy E = — 2 Jl cos(fc). We 
represent the wavefunction on the left leads as ^l(x) = 
= (e~ lKx 8L + e lKx h) and similarly the wavefunc- 



iKol 



yjl sin(iT) 

tion on the right leads as ^r(x) = — ^= 

e %Kx d). Here the components of the amplitude vector a 
correspond to Floquet channels, where mth entry is the 
amplitude of states with energy E + mQ. The matrix K 
is the the diagonal matrix whose diagonal vector gives 
the wavenumber k n = cos~ 1 ((E + nQ)/(— 2Jl)) f° r each 



Floquet state. The normalization of the state y 2 sin (if) 

ensures the unitarity of scattering matrix, i.e. SS^ = 1, 
defined below. Because the leads have finite bandwidth, 
the transmission of particles has a hard cutoff at a fi- 
nite Floquet level \n\ « 2Jl/Q. Throughout the pa- 
per, we take Jl to be much larger than J such that all 
the Floquet states occupied in the central system can be 
smoothly transmitted to the leads. Even when Jl is not 
much larger than J, this gives a good approximation if 
the transmission amplitudes t(E m ,E) for E m >2Jl are 
small. 

Then the amplitudes a, b, c, d are related through scat- 
tering matrix S = 1 ' ^ 



t r' 



as 

r t' 
i f' 



(6) 
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FIG. 5: The scattering process for the periodically driven 
systems. Incoming wave amplitudes and outgoing wave am- 
plitudes are represented as a vector whose rath entry is the 
amplitude of rath Floquet channel with energy E-\-mQ. Here 
we apply light to the system with size L. Amplitudes on the 
left lead is related to those on the right through the transfer 
matrix M, see Eq. This transfer matrix is directly related 
to the conductance. 



Computationally convenient, closely connected matrix is 
the transfer matrix M which relates the amplitudes on 
the left to the amplitudes on the right; 



M 



(?) 



It can be easily checked that M can be expressed through 
the elements of scattering matrix S as 



M = 



rt 



t' - rt 



(8) 



In the following, we relate the Hamiltonian to M, and 
obtain the expression of t through the relation above. 

It is straightforward to write down the local transfer 
matrix Mi which "transfer" the amplitudes of wavefunc- 
tion from site i to site z + 1. At each site z, the wavefunc- 
tion has components with different Floquet levels. The 
amplitude associated with each level is just the projec- 
tion of \ipE(m)} on site z, i.e. s m i = (z^^ra)). We 
represent the amplitude of wavefunction at site z as a 
vector whose mth component is s m i. If we write the 
Schrodinger equation Eq. ^ in terms of the amplitudes 
of wavefunction at each site s 7 -, then 



(E + n-Wi)si 



-J 1 s,_ 



+1 



in 



on 



in 



\ 
( 



j = 



Jo J\ J2 
J-l Jo J\ 
J-2 J-l Jo 



Here, J m = JB m (eAa) is a hopping with m photon ab- 
sorptions where B m (x) is mth Bessel function. Formally, 
this equation can be written in terms of the transfer ma- 
trix Mj as 



s*+i 
s,- 



Mi 



Si 

s*-i 



Mi = 



-j-^e + q- 
1 



Wi) 



-J- 1 J T 





If we suppose that the size of the central system to be 
L where disorder and driving fields exist only from site 
1 to site L, then the transfer matrix from sites (0, 1) to 
sites (L, L + l) is given by M 1L = M L ]jf~^ MiM x . Here 
the transfer matrix that connects the central system and 
left and right leads, M\ and Ml respectively, are given 
by 



Mi 
Ml 



-J-^E + Cl-Wi) -J~ x Jl 

1 

-J^iE + Cl-Wi) -Jl X J T 

1 



Note that whereas Mil connects the amplitudes at 
sites i=l, to sites \—L + 1, L, the transfer matrix M 
defined above connects amplitudes a, b,c,d of the free 
propagating mode in the leads from the left lead to the 
right lead. These two different amplitudes are related 
through the equations 



sl+i 



Q 




'2sm(K) 



e iK e ~iK 



1 



(9) 
(10) 

(11) 



Thus, the transfer matrix M is obtained through Mil as 

/ e -iKL q 



M 







iKL 



Q- X M 1L Q 



(12) 



-iKL 




The unitary transformation r>r has no ef- 

y e l 

feet on the physical transport properties, so in the fol- 
lowing, we study the equivalent transfer matrix M = 

In principle, it is now straightforward to calculate M 
to obtain the conductance of the driven system G = 
Go (iH) QQ . Numerically, this procedure faces difficulty 

because calculation of M gives the inverse of i and one 
has to take the inversion. Since states of disordered one 
dimensional system are localized, many transmission am- 
plitudes are exponentially small as a function of system 
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size. In order to have a reasonable precision in the com- 
putations of transmission amplitudes, it is necessary to 
keep the small numbers in the inverse matrix t~ x which 
gives a large contribution to the transmission amplitude 
(tH) QQ . However, since the transmission amplitudes are 
very small in the disordered regimes, it is generically nec- 
essary to maintain accuracies that retain exponentially 
small numbers. This faces a difficulty as one goes into 
deeper localized regimes. 

This problem arises since M is the multiplication of 
many matrices with disorder at each site, which cre- 
ates eigenvalues that are separated by exponentially large 
numbers. An ingenious method to avoid the prob- 
lem is to take the inverse before the multiplication^ . 
We aim to obtain the inverse of the matrix t~ x = 

(10) Q^MilQ ( J ). To this end, we first start 



from 



Now we rewrite 



di 



d\. In the next step, we multiply the transfer 



matrix Mi to obtain 



^ u±d 1 1 

t )=*.("f 

this procedure to obtain the equality 

(1 ^)Q- X M 1l q{^ =u L+1 d L ---d 1 
The transmission amplitude is then given by 



We repeat 



(13) 



t = d 1 1 d 2 



• dr 1 U 



(14) 



From the expression of Eq. (5j, the conductance and 
Floquet transmission amplitude are related by the 
equation 



G = G (tH) 



oo 



(15) 



where Go is the unit of conductance e 2 /h in the case 
of spinless electron in our study. As we have noted in 
the summary section, it is crucial that the the equation 
Eq. (15) is only 00 Floquet component of the Floquet 
transmission probability tH. This is in stark contrast 
with the static, multichannel case, where the conduc- 
tance is given by the trace of multi-channel transmission 
probability because incoming electrons can be in any of 
the channels present. In the Floquet channel case, incom- 
ing electrons are all sitting near the Fermi surface, which 
restrict the conductance to be 00 Floquet component of 

m. 



IV. RESULT 



A. Analysis of moments 



In Section [II C| the results indicate that the applica- 
tion of light not only changes the average conductance, 
but also significantly changes the shapes of conductance 



distributions. One way to characterize the shapes of dis- 
tributions is through their moments. Each moment cap- 
tures a particular aspect of the distributions, and we will 
analyze each moment in this section. 

The first moment (G/Go) gives the average conduc- 
tance, which we have argued, results from the competi- 
tion between dimensional increase and dynamic localiza- 
tion. Because of this competition, behaviors upon the in- 
crease of intensity of light are different from the increase 
of physical channels in the static systems. In Fig. [6] a) 
and b), we plot the color plot of average conductance in 
the log scale, Fig. [6] a) is the average conductance for 
static quasi- ID systems for various disorder strength and 
the number of physical channels. As is expected, the av- 
erage conductance increases as the number of channel in- 
creases and disorder strength decreases. There is a clear 
crossover from insulating to conducting behaviors in the 
diagonal of the color plot. This behavior should be con- 
trasted with Fig. |6]b), where the average conductance 
is plotted for various intensity of light A and disorder 
strength. In contrast with the static counterpart, the av- 
erage conductance does not show large increase, except 
for the strongly disordered regimes. For large intensity 
of light, the decrease of conductance due to dynamic lo- 
calization is recognizable. 

Second moments gives the fluctuation, given by 
((G/Gq - (G/G )) 2 ). As is evident from the distribu- 
tions plotted in Fig. [3j the application of light signifi- 
cantly reduces the fluctuation in the weakly disordered 
regimes. Motivated by this observation, we plot the fluc- 
tuations as a color plot in Fig. [6] c) and d) where color 
is assigned in the log scale. For comparison, we plot in 
Fig. [6]c) the strength of fluctuations of static quasi- ID 
systems for various disorder strength and the number of 
channels. Fluctuation generally increases as the number 
of channels increases because the average conductance 
also increases. In Fig. [6]d), we plot the strength of fluc- 
tuations of systems under the application of light for vari- 
ous intensity of light and disorder strength. In contrast to 
static systems, while the average conductance increases 
as the intensity of light is increased (see Fig.[3|, the fluc- 
tuation is decreased as the light intensity is increased 
from zero. 

Third moment leads to the skewness, which is the mea- 
sure of how symmetric the distribution is. Concretely, 
skewness is defined as 



{(G/G - (G/G )f ) 
{(G/G - (G/Go)) 2 ) 



2\3/2 



(16) 



In Fig.[6]e) f), we plot the skewness of the conductance in 
the log scale color plot for static and driven systems, re- 
spectively. Skewness can take either positive or negative 
value, and negative skewness means that the distribution 
has long tails on the left, and positive skewness means the 
distribution has long tails on the right. Here we plot the 
absolute value of the skewness in the log scale. 

In the static quasi-lD systems, skewness increases as 
the disorder strength is increased and it decreases as the 
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FIG. 6: Color plot of average, fluctuations and skewness (see Eq. (16))) for static quasi-lD systems and driven systems. Here, 
the color is assigned in the log scale. For static quasi- ID systems (top), the behavior of all moments is governed by the crossover 
from conducting regime (normal distribution with large average conductance) to insulating regime (log-normal distribution with 
small average conductance). Such crossover region in the diagonal line is apparent. On the other hand, driven systems display 
richer structure (bottom). The average conductance shows little change in the weakly disorder regimes as A is increased, it 
increases for small A and decreases for larger A in the strongly disordered regime. The fluctuation is strongly suppressed as A 
is increased from zero for weakly and intermediate disordered regimes. Skewness displays a similar crossover behavior as that 
of static quasi-lD systems, as is apparent also in Fig. [3] 



number of channels is increased, see Fig.[6]e). Both ef- 
fects can be understood as the result of transition from 
normal to log-normal distributions as the system goes 
from conducting regime to localized regime. The skew- 
ness of driven systems has the tendency to decrease as the 
intensity of light is increased and the disorder strength is 
decreased. The dependence on the disorder strength can 
be understood in a similar fashion as the static systems 
as the crossover between insulating regime and conduct- 
ing regime. On the other hand, the decrease of skewness 
on the increase of intensity has the larger effects than 
what is implied by the change of average conductance 
(also compare with Fig. |3|. 

For both the fluctuations and skewness, the behaviors 
of moments for the static quasi- ID systems are tied with 
the average conductance; the crossover behaviors from 
conducting regimes (left top) to the insulating regimes 



(right bottom) are apparent in any of three color plots of 
average, fluctuation and skewness. On the other hand, 
one of the main observations of our study for disordered 
systems under the application of light is the behaviors 
of average, fluctuation and skewness are not necessarily 
correlated; fluctuations and skewness are strongly sup- 
pressed by the application of light when the average con- 
ductance decreases or barely changes in the weak and 
intermediate disorder strength. 

B. Control of time- reversal symmetry through the 
driving shape 

In the previous sections, we consider the application of 
a simple harmonic drive. Such drive which possesses the 
time r such that theHamiltonian is symmetric around 
that point H(r — t)* = H(r + t) for all t can be shown 
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3,) time reversal symmetric drive 



b) time reversal symmetry 
breaking drive 




cos Qt + cos 2Qt cos Qt + cos(20 + 7r/2) 

FIG. 7: The shape of driving field can modify the symmetry 
class of the driven systems. Here, the driving of the form 
cosQt + cos (2Qt + ip) is considered, a) ip = leads to sym- 
metric shape of the drive, which implies the presence of time- 
reversal symmetry, b) when ip ^ or 7r, the drive breaks 
time-reversal symmetry. Here we plot the shape of the drive 
for (p = 7r/2. 



to be "time-reversal symmetric" when the dynamics is 
"averaged" over one period of time 34 . Here H* is the 



complex conjugation of the Hamiltonian and we assumed 
that the electrons are spin- less. Such time-reversal sym- 
metry plays an important role in the static disordered 
systems, and the presence and absence of time reversal 
symmetry changes the average as well as fluctuations of 
conductance, as we have described in Section!!! Al In this 



section, we discuss the possibility of changing the sym- 
metry class of the driven disordered systems by choosing 
the shape of the driving. 

It is possible to break the time reversal symmetry of 
the Hamiltonian by changing the shape of the driving, 
for example, by driving the systems with frequency Q 
and 2Vt with some phase difference ip. When ip = 0, the 
composite drive again possesses a symmetric time r, but 
for generic values of cp, there is no symmetric time r, see 
Fig. [7[ Intuitively, this is an indication of the lack of 
extra symmetry in the driving. 

From random matrix point of view, the difference be- 
tween the two drives can be seen as the difference in the 
symmetry class of Floquet Hamiltonian. First we rewrite 
Eq. (§ as a matrix equation for each Floquet level as 
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This equation can be thought of as the eigenvalue prob- 
lem for Floquet states and we refer to the matrix ap- 
pearing in the right hand side of Eq. (17) which mul- 



tiplies Floquet states as Floquet Hamiltonian Hp We 
consider the time reversal symmetry of the driven sys- 
tem as the straightforward extensions of time reversal 
symmetry of static systems applied to H^; time-reversal 
symmetry is the existence of anti-unitary operator T such 
that THfT' 1 = H F . 

If we represent the driving by light as the time- varying 
electric potential as in Appendix, then it is clear that the 
simple drive E cos(fit) leads to H F that is real, and there- 
fore, the system is time reversal symmetric with T = 1. 
Now we add the second driving with doubled frequency 
E cos(2Vtt + ip). When ip = 0, H F is again real. On the 
other hand, for ip =^ or 7r, Hf becomes complex due 
to the relative phase between Q drive and 2Q drive, and 
the time reversal symmetry is broken. Because the sym- 
metry class of Floquet Hamiltonian is different for these 
drives, the conductance property of these two drives are 
expected to be different as in the case of static systems. 

In Fig. [81 we observe the variation of average conduc- 



tance for disordered systems in the presence of two driv- 
ing frequencies Q and 2^ with various phase difference 
between the two. In terms of electric potential, the form 
of the driving is E {cosQt + cos (2Qt + <p)}. Here we 
took the intensity of light A = 0.3, Q = 0.01 J, W = 0.6 J, 
L = 601. In the calculation, we retained the number 
of Floquet states ±50 to ensure the convergence of the 
conductance. 6 can be thought of as "magnetic field" 
that breaks the time-reversal symmetry. Similar to the 
weak localization in static disordered systems, the driven 
systems displays negative "magnetoresistance," in which 
the breaking of time-reversal symmetry by non-zero 
increases the conductance in the system. 

On the other hand, the breaking of time reversal sym- 
metry seems to have little effects on the shape of the 
distribution. In Fig. [9j we give the conductance distribu- 
tions for disordered systems for the same parameter as 
in Fig. H for = and = tt/2. The distribution of the 
conductance is always Gaussian, irrespective of the value 
of 0. 

In Fig. [9J we give the conductance distributions for 
disordered systems with Q and 2^ drives with phase dif- 
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FIG. 8: The average conductance of disordered systems in 
the presence of two driving E {cos(Qt) + cos(2Qt + ip)} for 
various values of 0. The dot is the numerical result and the 
line is the guide to the eye. Here, L = 601, W — 0.6 J, 
£1 = 0.01J and intensity of light A = 0.3. Similar to the 
negative magnetoresistance of weak localization, the driven 
system also displays the increase of conductance under the 
breaking of time-reversal symmetry. 
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FIG. 9: The conductance distributions of disordered systems 
in the presence of the driving E {cos(Qt) + cos(2Qt + ip)} 
with ip = and ip = tt/2. Here, L = 601, W 0.6J, 
Q, = 0.0 1J and intensity of light A = 0.3. While the average 
value is increased by the breaking of time-reversal symmetry, 
the shape of the distributions remain Gaussian. 



ference cp = and cp = tt/2. Here we choose L = 601, 
W = 0.6J, Q = 0.01 J and intensity of light A = 0.3. To 
ensure the convergence of the conductance, we take up to 
±50 Floquet levels in the numerical computations. Con- 
ductance distributions of time-reversal symmetric drive 
and time-reversal symmetry breaking drive have almost 
identical shape, but their average values are different by 
about 20%, where time-reversal symmetry breaking drive 



leads to lower conductance. The difference in the aver- 
age conductance decreases as the intensity of drive A is 
decreased. 

While in one dimensional disordered systems we study 
here, the time-reversal symmetry breaking only led to a 
small change of conductance, the idea is general and can 
be applied to higher dimensions. For example, applica- 
tion of linearly polarized light to two dimensional system 
does not break time-reversal symmetry, but that of circu- 
larly polarized light breaks the time reversal Symmetrel 
Thus such control of symmetry class of driven systems is 
a promising approach to control the conductance prop- 
erties. 



C. Conductance dependence of a sample on the 
driving frequency and intensity: 
giant opto-response of conductance. 

1. Frequency dependence 

In the previous sections, we analyze the conductance 
distribution of disordered systems under the application 
of light. In experiments, it might not be easy to create 
so many samples with different realization of disorders, 
and question arises as to how one can measure such con- 
ductance distributions. In the static, weakly disordered 
systems, one aspect of conductance distribution, n amely 
fluctuation, is probed as a function of magnetic fiel d 21 1 35 1 . 
While such study is, in principle, different from the study 
of fluctuations for different realizations of disorders, it 
probes inherent fluctuation structure built in the eigen- 
energy distributions of disordered systems. As a result, 
the study of fluctuations for both different realizations of 
disorders and intensity of magnetic field displays univer- 
sal conductance fluctuation J 2 1 1 22 1 

In disordered systems in dimensions higher than one, 
the application of magnetic field might be used to study 
fluctuations. However, in one dimensional system, there 
is no orbital effect of magnetic field on the conductance. 
If we suppose that the wire is running along x direction 
and magnetic field is pointing in z direction, the opera- 
tion y — >• — y reverses the direction of magnetic field, but 
does not change the current for spin- less electrons, and 
thus we see that magnetic field does not have any effect 
on conductance. 

Here we propose the study of conductance fluctuations 
of disordered systems under the application of light as a 
function of driving frequency. We note that such study is 
essentially different from the study of conductance fluc- 
tuations for different realizations of disorders, but yet 
probes inherent fluctuation structure in the disordered 
systems. Apart from the fundamental question of con- 
ductance fluctuations, this study is also interesting from 
application point of view as a knob to change the con- 
ductance of disordered systems. 

In strongly disordered systems under the application 
of light, localized electrons can hop to neighboring local- 
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FIG. 10: Illustration of the photon-assisted hopping in the 
strongly disordered systems. Each state with a given energy 
E is spatially localized. If two states are spatially overlapping, 
and moreover two states are roughly separated by Q where Q 
is the driving frequency, the electrons can hop from one state 
to the other. The energy of each state is widened by the in- 
tensity of light J A, where A is the intensity of light and J is 
the static hopping amplitude. Changing of driving frequen- 
cies from, say, uji to CJ2 in this picture, leads to coupling of 
different states, resulting in the large change of conductance 
for a given realization of disorders. 



ized states with the help of photon energy, see Fig. [To] In 
this picture, the accessible states are determined by two 
conditions; 1. the spatial overlap of localized states; 2. 
energy separation between the two localized states. For 
the second condition, two states need to be separated 
by the frequency of light Q. This resonant condition is 
widened due to the intensity of light, where the widening 
energy is given by Rabi frequency of drive roughly given 
by J A for weak intensity of light. Thus the electron can 
hop from one state to the other if two states are sepa- 
rated by Q — J A ~ ft + J A. While the change of the 
frequency of drive by a small amount barely changes the 
conductance averaged over different realizations of disor- 
ders, such change can dramatically modify the accessible 
states for a given realization of disorder. 



In Fig. 11, we plot a dependence of conductance for 
a given realization of disorder for various frequency of 
drive. Here we chose L = 601 with two channels as be- 
fore, with intermediate disorder strength W = 0.6 J and 
small intensity of light A = 0.02. Despite of small in- 
tensity of light, the conductance strongly fluctuate as a 
function of driving frequency Q. The fluctuation is ex- 
pected to occur at the typical energy spacing of the sys- 
tem, which can be roughly estimated as AE = J/L w 
1.6 x 10 -3 J. Such giant response of conductance on the 
small change of driving frequency could provide an inter- 
esting device. Note we expect that that the amplitude 
of fluctuation as a function of driving frequency goes to 




0.01 0.015 0.02 0.025 0.03 
Driving Frequency (0/ J) 

FIG. 11: Dependence of conductance on driving frequency 
for the system with a given realization of disorder. Here we 
chose L = 601 with two channels as before, with intermediate 
disorder strength W = 0.6 J and small intensity of light A = 
0.02. We observe a giant response of the conductance due to 
the change of frequency by a small amount. 



zero as the intensity of light A goes to zero. Thus for 
very small intensity of light, the fluctuation amplitude 
will be suppressed. 



2. Conductance dependence for varying intensity of light 



In this paper, we study the dependence of conductance 
distributions on the intensity of light as a general charac- 
terization of statistical ensemble of disordered systems. 
From the point of view of controlling the conductance 
of disordered systems, yet another interesting question is 
the dependence of a sample with a given disorder on the 
intensity of light. As we have explained above, change of 
the intensity of light can alter which and how each local- 
ized states are coupled with each other. This, in return, 
can change the conductance. 

In Fig. [l2j we plot the behavior of conductance for a 
given realization of disorder as a function of the intensity 
of light, A. We choose L = 601 with two channels with 
intermediate disorder strength W = 0.6 J and small driv- 
ing frequency of light Q = 0.01 J. Since the intensity of 
light A controls the width of resonant condition, the ef- 
fect is expected to be more smooth than that of frequency 



change. Yet, we see in Fig. [12] that the conductance does 
change by a large amount for a small change of A. Thus 
this parameter can also be used to induce the giant op- 
tical change of conductance for disordered systems. 
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0.01 0.02 0.03 0.04 0.05 
j\ (Intensity of light) 

FIG. 12: Dependence of conductance on the intensity of light, 
A for the system with a given realization of disorder. We 
choose L = 601 with two channels with intermediate disorder 
strength W = 0.6 J and small driving frequency of light Q = 
0.01J. Conductance has a smooth dependence on A but yet, 
we observe a large change of conductance for a small change 
of A 



V. CONCLUSION 

In this paper, we investigated the behavior of full con- 
ductance distributions of disordered quasi- ID systems 
under the application of light. By developing the Flo- 
quet transfer matrix method, we study the effect of ap- 
plication of low frequency light for large systems, where 
photon-assisted transport plays an important role in the 
conductance. The comprehensive study of average, fluc- 
tuations, skewness and shape of conductance distribu- 
tions showed that the application of light influences the 
conductance in different ways for weakly, intermediate 
and strongly disordered regimes. We argued that the 
conductance is determined by the competition between 
dynamic localization which dominates in the weakly lo- 
calized regime and dimensional increase which dominates 
in the strongly localized regimes. Here we confirmed the 
previous result!^ which showed the elongation of localiza- 
tion length upon the application of light, and saturation 
of such increase for larger intensity of light in the strongly 
localized regime. In the intermediate disorder strength, 
the application of light results in the crossover behavior 
from log-normal like distribution to normal distribution 
of conductance. 

Many of the ideas demonstrated in this work can be 
directly extended to higher dimensional disordered sys- 
tems. The physics of localization is known to strongly 
depend on the dimensionality of the systems, where in 
1 and 2D, all the states are localized and in 3D, it is 
believed that the states are localized below a mobility 



edge and extended above it. Thus the effect of dimen- 
sional increase is expected to have a large effect in two 
dimensional systems under the application of light. This 
physics will be studied in the future work. 

In addition, it is interesting to study the disordered two 
dimensional systems in different geometrical structure. 
For example, hexagonal lattice structure has Dirac cones 
at the half filling, and intriguing weak anti-localization 
behavior has been observed in graphene 5 . Linearly po- 
larized light does not destroy the Dirac cones, and it is 
of interests to see how such anti-localization is modified 
under the application of light. 

Furthermore, the application of circularly polarized 
light on hexagonal lattice breaks time-reversal symmetry, 
and is predicted to give rise to band structures with non- 
trivial topology which has non-zero Chern numbers^. 
Such non-zero Chern number is known to be the origin of 
integer quantum Hall effect and chiral edge states. Thus 
in the disordered hexagonal lattice systems under the ap- 
plication of light, one can study the rich physics of the 
interplay between the localization and topological band 
structures. 
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Halperin, Jay D Sau, Gil Rafael, Chris R. Laumann, and 
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Army Research Office with funding from the DARPA 
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DMR-07-05472, AFOSR Quantum Simulation MURI, 
the ARO-MURI on Atomtronics. 



Appendix A: Alternative gauge: 
sinusoidally varying electric field 

As we have mentioned in the main text, it is sometimes 
more convenient to use the equivalent description of the 
light application through a different gauge choice. If we 
represent the light field as AC electric field, then the 
Hamiltonian takes the form 

<ij> i i 

(Al) 

In this gauge, Hamiltonian has only the first harmonics of 
driving frequency Q with H(l) = H( — l) = ^ Y^i x i4 c i- 
In a finite system, this operator is well-defined, but due 
to the presence of potential proportional to the inter- 
pretation of the operator becomes subtle in the thermo- 
dynamic limit. Therefore, this gauge is computationally 
useful in a finite system, but the other gauge choice we 
have used in Eq. ([!]) has clearer conceptual understand- 
ings. 
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In this gauge, the Schrodinger equation takes the form 



(E + Q-hi)si = -J(si_i +s i+ i) 
/ • 



10 



on 



-in 



hi = 



Wi Exi/2 
Exi/2 Wi Exi/2 
Exi/2 Wi 



V 



leading to the local transfer matrix 

= Mi 



Si+l 
Si 



Si 
Si-1 



Mi 



-(E + Q-hi)/J -1 



1 







The transfer matrices that connects the left lead to the 
central system Mi and the central system to the right 
lead Ml are given by 



Mi = 
Ml = 



-(E + n-hi)/j 

1 J 

-(E + Cl-hi)/JL -J/Jl 

1 



As in the other gauge, the transfer matrix M is now 
given by M = Q _1 Y\ i MiQ up to the unitary rotation. 
The calculation of the transmission matrix t will be anal- 
ogous as describedi n the main text. 
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